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TITLE 

DATA ENCODING WITH AN AMPLITUDE MODEL AND PATH 
BETWEEN THE DATA AND CORRESPONDING DECODING 

BACKGROUND OF THE INVENTION 

Field of the Invention 

[0001] The present invention generally relates to digital signal encoding and for 
this purpose provides a method and a device for encoding a digital signal. 
[0002] The object of the encoding is to compress the signal, which makes it 
possible respectively to transmit the digital signal and to store it in memory while 
reducing the transmission time, or the transmission rate respectively, by reducing 
the space in memory that is used. 

[0003] The invention is in the field of lossy compression of digital signals. The 
digital signals considered here are of any type, for example fixed images, video, 
sound, or computer data. 

Related Art 

[0004] In what follows, the encoding and decoding of a fixed image will more 
particularly be considered. 



-2- 



[0005] In this context, certain encoding modes use a path established between a set 
of digital samples. For example, French patent applications No. 01 06933, 01 
12064 and 01 13922 concern such encoding modes. 

[0006] For the encoding to be efficient, that is to say for it to have a good 
rate-distortion ratio, it is necessary to determine the path in an appropriate manner. 
[0007] Techniques exist to determine a path between a set of samples. These 
techniques are known as techniques for solving the traveling salesman problem. A 
review of these techniques is for example set out in the work of Gerhard Reinelt 
entitled "The traveling salesman, computational solutions for TSP applications", 
Springer- Verlag, 1994. 

[0008] Generally it is desired that the result of the encoding be as compact as 
possible. Thus, in the case of encoding using a path established between a set of 
digital samples, the path is described by the shortest possible path (traveling 
salesman problem). The shorter the displacement vectors linking each point of the 
path, the more compact will be the result of their encoding. 



SUMMARY OF THE INVENTION 
[0009] In a first aspect, the present invention provides a method and a device for 
determining a path between samples of a signal which make it possible to obtain a 
description of that path which is still more compact than with the known methods. 
[0010] To that end, the invention provides a method of encoding a set of data 
representing physical quantities, the encoding comprising forming a path between 
the samples of the set, the path passing a maximum of once via each sample 
location, characterized in that it comprises, for two given successive sample 
locations of the path, the step of: 

forming a vector between the two sample locations, only taking into 
account the intermediate sample locations not already encoded by a vector. 
[0011] The invention enables a description of that path to be obtained which is 
even more compact than with the known methods. This is because the vectors 
formed according to the invention are shorter than according to the prior art. 
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[0012] According to a preferred feature, the set of data being bidimensional, the 
method comprises the prior transformation of the set of data into a second set of 
monodimensional data. 

[0013] According to an alternative preferred feature, the set of data being 
bidimensional, a vector is decomposed into two vectors and one of them is formed 
only taking into account the intermediate sample locations not already encoded by 
a vector. 

[0014] According to a preferred feature, forming a vector between the two sample 
locations comprises: 

testing all the intermediate locations between the two sample locations to 
determine whether they have already been encoded, 

reducing the length of the vector as a function of the number of locations 
already encoded. 

[0015] Thus, the information already encoded is exploited in order to reduce the 
quantity of information necessary to form a path. 

[0016] According to a preferred feature, the set of data is a block formed in an 
initial set of data. 

[0017] According to a preferred feature, the set of data is the result of a DCT 
transformation of an initial set of data. 

[0018] The invention also relates to a method of decoding a set of data encoded by 
the method presented above, comprising a step of decoding the path, characterized 
in that the decoding of the path comprises the step of: 

deducing the location of a coefficient of the path as a function of the 
preceding location, the vector between the two locations, and the locations already 
decoded. 

[0019] In a complementary manner, the invention relates to a device for encoding 
a set of data representing physical quantities, comprising means for forming a path 
between the samples of the set, the path passing a maximum of once via each 
location, characterized in that it comprises: 
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means for forming a vector between two given successive locations of the 
path, only taking into account the intermediate samples not already encoded by a 
vector. 

[0020] The encoding device comprising means for implementing the encoding 
features presented above. 

[0021] The invention also concerns a decoding device comprising means for 
implementing the decoding features presented above. 

[0022] The encoding device and the decoding method and device have similar 
advantages to those already presented. 

[0023] The techniques for solving the traveling salesman problem use 
evolutionary calculation, which provides an optimum solution to the problem but 
requires a long calculation time. 

[0024] In a second aspect, the present invention aims to remedy the drawbacks of 
the prior art, by providing a method and a device for encoding digital data by 
which a path is determined more rapidly than with the prior art. 
[0025] To that end, the invention provides a method of encoding digital samples 
of a set of data representing physical quantities, the encoding including the 
determination of an amplitude model and of a path between the samples of the set, 
characterized in that it comprises the steps of: 

determining a number of samples to encode, 

constructing a list comprising the determined number of samples, classified 
by decreasing amplitude. 

[0026] The invention enables the path between the samples to be determined more 
rapidly than according to the prior art. The deterioration of the compression 
performances (rate and distortion) is low and generally negligible. 
[0027] According to a first preferred variant, the method comprises the steps of: 
determining an initial list of samples, 

calculating an encoding cost as a function of the list of samples, 
modifying the list of samples, 

the steps of calculating and modifying being reiterated to find a minimum 
encoding cost. 
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[0028] Thus, according to this variant, a list of samples is determined. 
[0029] According to a preferred feature, the method further comprises the step of 
encoding the set of data on the basis of the list of samples which provides the 
minimum encoding cost. 

[0030] The list determined according to the invention is thus used to encode the 
data. 

[0031] According to a preferred feature, the initial list of samples comprises all the 
samples of the set of data. 

[0032] According to a preferred feature, the modification of the list of samples 
comprises the withdrawal of the sample of least amplitude. 

[0033] Thus, the lists envisaged will have a number of samples that decreases as 
the iterations progress. 

[0034] According to preferred features which may be combined, the encoding cost 
comprises the rate of the encoded data and the encoding cost comprises the 
distortion of the encoded data. The rate and distortion are commonly taken into 
account on compression of data. Among the lists envisaged the one which has the 
least encoding cost is selected to encode the data. The determination of that list is 
rapid, and its construction gives it low encoding cost. 

[0035] According to a second variant, the invention relates to a method as 
presented earlier, comprising an initialization of an evolutionary algorithm 
according to which a population of lists of samples is determined, the population 
comprising a predetermined number of lists, characterized in that the determination 
of the population comprises the steps of: 

determining a first list of samples classified by decreasing amplitude, 
modifying the first list by withdrawal of a predetermined number of 
samples of lowest amplitude, to form a second list, 

the steps of determining and modifying being reiterated by taking the 
second list of an iteration as the first list for the following iteration, provided that 
the predetermined number of lists has not been reached and that the second list has 
a non-zero number of samples. 
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[0036] According to this variant, it is the initialization of the evolutionary 
calculation which is simplified and rendered more rapid to execute. Furthermore, 
the lists created according to the invention are generally better than lists picked 
randomly, which means that convergence of the evolutionary calculation is 
favored. 

[0037] According to a preferred feature, the population is completed by lists 
picked randomly, if the second list formed has a zero number of samples before the 
predetermined number of lists has been reached. 

[0038] According to a preferred feature, the set of data is a block of samples 
formed in a larger set of data. 

[0039] According to a preferred feature, the data are a digital image. 
[0040] In a complementary manner the invention relates to a device for encoding 
digital samples of a set of data representing physical quantities, comprising means 
for determining an amplitude model and a path between the samples of the set, 
characterized in that it comprises: 

means for determining a number of samples to encode, 

means for constructing a list comprising the determined number of samples, 
classified by decreasing amplitude. 

[0041] The device according to the invention comprises means for implementing 
the features presented above and has similar advantages to those presented above. 
[0042] Furthermore, it is desirable for the calculations to determine the amplitude 
model to be neither too long nor too complex. 

[0043] In a third aspect, the object of the present invention is to provide a method 
and device according to which the encoding of the data by an amplitude model and 
a path between the data is simplified. 

[0044] To that end, the invention provides a method of encoding a set of data 
representing physical quantities, the set of data comprising coefficients, the method 
comprising determining an amplitude model of the coefficients and a path between 
the coefficients, characterized in that it comprises the prior step of: 

putting the coefficients in order as a function of their respective locations in 
the set of data. 
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[0045] The invention makes it possible to encode the data using an amplitude 
model and a path formed between the data more rapidly than according to the prior 
art. 

[0046] According to a preferred feature, the prior step of putting the coefficients in 
order comprises putting them in order following a zig-zag path. This type of path 
is simple to implement. 

[0047] According to a preferred feature, the ordered coefficients are arranged in a 
one-dimensional vector. The path between the coefficients will then be formed in 
that vector. 

[0048] According to an alternative feature, the prior step of putting the coefficients 
in order comprises forming blocks of predetermined size in the set of data and 
putting the blocks formed in order according to a predetermined order. 
[0049] According to a preferred feature, the ordered coefficients are arranged in a 
three-dimensional table. According to this variant, the path between the 
coefficients is then formed in this table. 

[0050] According to a preferred feature, determining the amplitude model 
comprises classifying the coefficients by decreasing amplitude. This classification 
is simplified by the prior putting in order of the coefficients according to the 
invention. 

[0051] According to a preferred feature, the set of data is the result of a discrete 
cosine transformation of an initial set of data. 

[0052] According to an alternative feature, the set of data is the result of a discrete 
wavelet transformation of an initial set of data. 

[0053] The advantage of these transformations is to concentrate the most 
significant coefficients in the same region of the signal. 
[0054] The invention also relates to a method of decoding a set of data 
representing physical quantities encoded by the method presented earlier, 
comprising the decoding of an amplitude model of the coefficients and of a path 
between the coefficients, so as to form a first set of decoded coefficients, 
characterized in that it comprises the step of: 
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putting the decoded coefficients in order as a function of their respective 
locations in the first set of decoded coefficients so as to form a second decoded set. 
[0055] According to a preferred feature, the first set of decoded coefficients is a 
one-dimensional vector and the decoded coefficients are taken from first to last and 
are arranged in a two-dimensional table following a zig-zag path. 
[0056] According to an alternative feature, the first set of decoded coefficients is a 
three-dimensional table, and the decoded coefficients are arranged in a two 
two-dimensional table, the levels of the three-dimensional table being arranged in 
lexicographical order. 

[0057] In a complementary manner, the invention relates to a device for encoding 
a set of data representing physical quantities, the set of data comprising 
coefficients, the device comprising means for determining an amplitude model of 
the coefficients and a path between the coefficients, characterized in that it 
comprises: 

prior means for putting the coefficients in order as a function of their 
respective locations in the set of data. 

[0058] The invention also relates to a device for decoding a set of data 
representing physical quantities encoded by the device presented earlier, 
comprising means for the decoding of an amplitude model of the coefficients and 
of a path between the coefficients, so as to form a first set of decoded coefficients, 
characterized in that it comprises: 

means for putting the decoded coefficients in order as a function of their 
respective locations in the first set of decoded coefficients so as to form a second 
decoded set. 

[0059] The devices according to the invention comprise means for implementing 
the features already set out. 

[0060] The encoding device and the decoding method and device have similar 
advantages to those already presented. 

[0061] The invention also relates to a digital apparatus including the device 
according to the invention or means for implementing the method according to the 
invention. This digital apparatus is for example a digital camera, a digital 
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camcorder, a scanner, a printer, a photocopier, or a fax machine. The advantages 
of the device and of the digital apparatus are identical to those already set out. 
[0062] The invention also relates to an information storage means, which can be 
read by a computer or microprocessor, integrated or not into the device, and 
possibly removable, storing a program implementing the method according to the 
invention. 

[0063] The invention also relates to a computer program readable by a 
microprocessor and comprising one or more sequences of instructions capable of 
implementing the methods according to the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0064] The characteristics and advantages of the present invention will appear 
more clearly from a reading of three preferred embodiments illustrated by the 
accompanying drawings, in which: 

[0065] Figure 1 is an embodiment of a device implementing the invention, 
[0066] Figure 2 represents an encoding device and a corresponding decoding 
device, according to a first embodiment of the invention, 

[0067] Figure 3 represents an embodiment of an encoding method according to the 
first embodiment of the invention, 

[0068] Figure 4 represents an amplitude model used according to the present 
invention, 

[0069] Figure 5 represents a first variant of encoding locations of coefficients 
implemented according to the first embodiment of the invention, 
[0070] Figure 6 represents an example of a signal to be processed according to the 
first embodiment of the invention, 

[0071] Figure 7 represents an example of a signal to be processed according to the 
first embodiment of the invention, 

[0072] Figure 8 represents an example of a path to encode according to the first 
embodiment of the invention, 

[0073] Figures 9 and 10 represent a vector reduction implemented according to the 
first embodiment of the invention, 
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[0074] Figure 1 1 represents a second variant of encoding locations of coefficients 
implemented according to the first embodiment of the invention, 
[0075] Figure 12 represents an example of a path to encode according to the first 
embodiment of the invention, 

[0076] Figures 13 and 14 represent a vector reduction implemented according to 
the first embodiment of the invention, 

[0077] Figure 15 represents an embodiment of a decoding method according to the 
first embodiment of the invention, 

[0078] Figure 16 represents a first variant of decoding locations of coefficients 

implemented according to the first embodiment of the invention, 

[0079] Figure 17 represents a second variant of decoding locations of coefficients 

implemented according to the first embodiment of the invention, 

[0080] Figure 18 represents an encoding device and a corresponding decoding 

device, according to a second embodiment of the invention, 

[0081] Figure 19 represents a first variant of an encoding method according to the 
second embodiment of the invention, 

[0082] Figure 20 represents a second variant of an encoding method according to 
the second embodiment of the invention, 

[0083] Figure 21 shows an encoding device and a corresponding decoding device, 
according to the third embodiment of the invention, 

[0084] Figure 22 shows an embodiment of an encoding method according to the 
third embodiment of the invention, 

[0085] Figure 23 represents a first variant for putting the coefficients in order as a 
function of their locations, 

[0086] Figure 24 represents a second variant for putting the coefficients in order as 
a function of their locations, 

[0087] Figure 25 represents an embodiment of a decoding method according to the 
third embodiment of the invention. 

[0088] Figure 26 represents a first variant for putting the coefficients in order as a 
function of their locations, 
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[0089] Figure 27 represents a second variant for putting the coefficients in order as 
a function of their locations. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0090] According to the chosen embodiment shown in Figure 1, a device 
implementing the invention is for example a microcomputer 10 connected to 
different peripherals, for example a digital moving picture camera 107 (or a 
scanner, or any means of acquiring or storing an image) connected to a graphics 
card and supplying information to be processed according to the invention. 
[0091] The device 10 comprises a communication interface 1 12 connected to a 
network 113 able to transmit digital data to be processed or conversely to transmit 
data processed by the device. The device 10 also comprises a storage means 108 
such as a hard disk. It also comprises a drive 109 for a disk 1 10. This disk 1 10 
may for example be a diskette, a CD-ROM, or a DVD-ROM. The disk 110 like the 
hard disk 108 can contain data processed according to the invention as well as the 
program or programs implementing the invention which, once read by the device 
10, will be stored on the hard disk 108. According to a variant, the program 
enabling the device to implement the invention can be stored in read-only memory 
102 (referred to as ROM in the drawing). In a second variant, the program can be 
received in order to be stored in an identical fashion to that described previously 
via the communication network 113. 

[0092] The device 10 is connected to a microphone 111. The data to be processed 
according to the invention will in this case be audio signal. 

[0093] This same device has a screen 104 for displaying the data to be processed 
or serving as an interface with the user, who can thus parameterize certain 
processing modes, using the keyboard 1 14 or any other means (a mouse for 
example). 

[0094] The central processing unit 100 (referred to as CPU in the drawing) 
executes the instructions relating to the implementation of the invention, which are 
stored in the read-only memory 102 or in the other storage means. On powering 
up, the processing programs stored in a non-volatile memory, for example the 
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ROM 102, are transferred into the random access memory RAM 103, which will 
then contain the executable code of the invention, as well as registers for storing 
the variables necessary for implementing the invention. 

[0095] In more general terms, an information storage means, which can be read by 
a computer or microprocessor, integrated or not into the device, and which may 
possibly be removable, stores a program implementing the method according to the 
invention. 

[0096] The communication bus 101 affords communication between the different 
elements included in the microcomputer 10 or connected to it. The representation 
of the bus 101 is not limiting and, in particular, the central processing unit 100 is 
able to communicate instructions to any component of the microcomputer 10 
directly or by means of another element of the microcomputer 10. 
[0097] Figures 2 to 17 concern the first embodiment of the invention. 
[0098] With reference to Figure 2, the first embodiment of an encoding device 
according to the invention is adapted to encode a digital signal with the object of 
compressing it. The encoding device is integrated into an apparatus, which is for 
example a digital camera, a digital camcorder, a scanner, a printer, a photocopier, a 
fax machine, a database management system, or a computer. 

[0099] An image source 1 supplies a digital image IM1 to the encoding device 2, 
of which the operation will be detailed below. The encoding includes determining 
a path between the samples of the set. 

[0100] The encoding device comprises a transformation circuit 21 connected to a 
circuit 22 for calculating an amplitude model and path. The transformation circuit 
21 supplies a transformed image to circuit 22. 

[0101] The circuit 22 is connected to an encoding circuit 23 which performs an 
entropy encoding of the previously determined path. 

[0102] The operation of the encoding device will be detailed hereinafter. The 
determination of a path between the data to encode, which the present invention 
relates to, will be more particularly detailed. 
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[0103] The encoding device provides a file containing data representing the 
compressed image to means 3 for transmission and/or storage. These means are 
conventional and will not be described here. 

[0104] The means 3 are connected to a decoding device 4 whose operation will be 
detailed below. 

[0105] The decoding device comprises a circuit 41 for entropy decoding of the 
data which it receives. 

[0106] Circuit 41 is connected to a circuit 42 for decoding the amplitude model 
and the path, which associates an amplitude with each decoded location. 
[0107] Circuit 42 is connected to a circuit 43 for inverse transformation which 
supplies a decoded image IM2 to a device 5 for viewing the image. 
[0108] It should be noted that the encoding device and decoding device may be 
incorporated into the same apparatus, for example the computer 10 of Figure 1. 
[0109] Figure 3 represents an embodiment of a method of encoding an image, 
according to the first embodiment of the invention. This method is implemented in 
the encoding device and comprises the steps SI to S7. 

[0110] The method generally comprises a transformation of the signal to be 
encoded, then the determination of an amplitude model of the coefficients output 
from the transformation. The locations of these coefficients are next encoded 
according to a method using a path established between the coefficients. 
[0111] Such an encoding method is described for example in French patent 
application No. 01 06933. 

[0112] The method is carried out in the form of an algorithm which can be stored 
in whole or in part in any means of information storage capable of cooperating with 
the microprocessor. This storage means is readable by a computer or by a 
microprocessor. The storage means is integrated or not into the device, and may be 
removable. For example, it may comprise a magnetic tape, a diskette or a 
CD-ROM (fixed memory compact disk). 

[0113] Step SI is a linear or non-linear transformation of a digital image IM1 to be 
processed according to the invention. 
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[0114] In the preferred embodiment of the invention, the image is divided into 
blocks and then a discrete cosine transformation (DCT) is applied to each block, 
such as that applied under the JPEG standard. 

[0115] In variant form, another transformation is used, for example a discrete 
wavelet transformation, as in the JPEG2000 standard. 

[0116] The following processing is performed block by block. In variant form, the 

processing may be applied globally to the whole image. 

[0117] The following step S2 is an initialization at which a first block is 

considered. 

[0118] The following step S3 consists of classifying the coefficients of the current 
block according to decreasing amplitude. This results in a list P of coefficients. 
[0119] The following step S4 is the determination of an amplitude model. For 
this, a function for approximating the series of the classified coefficients is 
determined. This function is for example a decreasing exponential defined by a set 
of parameters which are determined by regression. French patent application 
No. 01 06933 describes this step in detail, as follows. 

[0120] The method will choose, from the series of P coefficients arranged from the 
greatest to the least and named by the series x l5 x 2 , x 3 , . . . , x p , a set A comprising a 
finite number of R parameters a l9 a 2 , . . . , B R . For example, a decreasing 
exponential family could be chosen, written f A (x)=a 1 .exp(a 3 .x+a 4 )a 2 . 
[0121] The method will then determine among a family of real functions that 
which reproduces as accurately as possible the evolution of the values of the Xj as a 
function of their rank i. 

[0122] At these steps, the parameters a l9 a 2 , a 3 , and a 4 will be determined such that 
the function f A is the closest possible to the points (1, x x ), (2, x 2 ), . . . ,(P, x p ). This 
typical problem is called function regression. Numerous techniques exist for 
finding the best set of parameters. For example minimization of the root mean 
square error RMSE could be sought, equal to 
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[0123] To achieve this, it could be chosen to quantize each parameter a t into a 
finite number of possible values. There are then a finite number of different values 
of parameters of A. The RMSE value is then calculated for each set of parameters 
and is then stored with the value of the parameters associated. When all the 
parameters have been incremented, the set which minimizes RMSE from among 
the stored parameters will be chosen. 

[0124] Figure 4 represents an example of an amplitude model A. A value A(k) 
provided by the amplitude model corresponds to each integer value k along the 
x-axis. The value A(k) is an approximation of the amplitude of the k* coefficient 
classified in decreasing order. 

[0125] The following step S5 is the encoding of the locations of the coefficients of 
the current block. 

[0126] Two variants of this step are detailed below. 

[0127] The following step S6 is a test to determine whether the current block is the 
last block of the image to encode. 

[0128] If the response is negative, that step is followed by the step S7 at which a 
following block is considered. Step S7 is followed by the previously described 
step S3. 

[0129] If the response is positive at step S6, then the encoding of the image is 
terminated. 

[0130] The first variant of step S5 of encoding the locations of the coefficients of 
the current block is detailed in Figure 5 in the form of an algorithm comprising 
steps S50 to S67. 

[0131] Step S50 is a transformation of the bidimensional signal into a 
monodimensional signal. For this, a zig-zag path is followed through the block B. 
As represented in Figure 6, in the block B of 8x8 coefficients, the diagonals of the 
coefficients oriented from bottom left to top right are considered. The path begins 
at the top left. The path passes through the diagonals either from bottom left to top 
right, or in the opposite sense. For a given diagonal, the direction of the path is the 
direction opposite to that of the preceding diagonal. 
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[0132] Following this path direction causes the coefficients to be arranged in a 
monodimensional vector V containing 64 coefficients. 

[0133] In variant form, as represented in Figure 7, block B is passed through line 
by line, from left to right and from top to bottom, in order to form a 
monodimensional signal. 

[0134] The following step S51 is the determination of a path in the vector V. 
[0135] A path is determined by an initial coefficient and the list of the vectors 
joining the other coefficients. Each coefficient of the path different from the initial 
coefficient is represented by a vector describing its location with respect to the 
preceding coefficient in the path. It should be noted that the path does not 
necessarily include all the coefficients of the current vector. This is because it is 
possible to encode only part of the coefficients and to set the other coefficients to 
the value zero at the time of later encoding. 

[0136] For example, Figure 8 represents a path comprising ten coefficients of 

which the locations x x to x 10 are connected by nine vectors V x to V 9 . 

[0137] In what follows, the locations {x 1? x 2 , . . ., x n } of coefficients in the vector 

V are considered, where n is an integer equal to the number of coefficients of the 

path. 

[0138] The following step S52 is an initialization at which a binary encoding 
parameter AC X is associated with each of the locations x of the vector V. All the 
encoding parameters AC X are set to the value zero at step S52. 
[0139] For the following processing, the encoding parameter AC X has the value 
zero if the location x has not yet been encoded, and has the value one once the 
location x has been encoded. 

[0140] The following step S53 is an encoding of the first location x 1 of the path. 
The parameter AC xl which is associated with it is set to the value one. 
[0141] The following step S54 is an initialization at which a parameter i is 
initialized to the value one. 

[0142] At the following step S55, the vector Vj = x i+1 -x i is considered. A location 
parameter x is set to the value x { . 
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[0143] The following step S56 is a test to determine if the current vector Vj is 
oriented in one direction or the other. For example an axis of coordinates of which 
the origin is at x x is considered. Along that axis, the vector V l (Figure 8) is 
positive whereas the vector V 2 is negative. 

[0144] When the current vector Y i is positive, step S56 is followed by step S57 at 
which the location parameter x is incremented by one unit, in order to consider the 
location immediately to the right of the current location. 
[0145] The following step S58 is a test to determine whether the encoding 
parameter of location x has the value one, that is to say whether location x has 
already been encoded. 

[0146] If the response is positive, then step S58 is followed by step S59 at which 
the value of the current vector V i is decremented by one unit, so as to "skip" the 
location already encoded. 

[0147] If the response is negative at step S58, this step is followed by the step S60. 
In the same way, step S59 is followed by the step S60. 

[0148] Step S60 is a test to determine whether the current value of the location 
parameter x is equal to x i+1 . 

[0149] If the response is negative, that means that all the locations between the 
positions x { and x i+1 have not yet been considered. Step S60 is then followed by the 
previously described step S57. 

[0150] If the response is positive at step S60, that means that all the locations 
between the positions Xj and x i+1 have been tested in order to know whether they 
have already been encoded or not. The length of vector V { has been reduced by the 
number of locations met that were already encoded. 

[0151] When the current vector Vj is negative, step S56 is followed by step S61 at 
which the location parameter x is decremented by one unit, in order to consider the 
location immediately to the left of the current location. 
[0152] The following step S62 is a test to determine whether the encoding 
parameter of location x has the value one, that is to say whether location x has 
already been encoded. 



- 18- 



[0153] If the response is positive, then step S62 is followed by step S63 at which 
the value of the current vector V l is incremented by one unit, so as to "skip" the 
location already encoded. 

[0154] If the response is negative at step S62, this step is followed by the step S64. 
In the same way, step S63 is followed by the step S64. 

[0155] Step S64 is a test to determine whether the current value of the location 
parameter x is equal to x i+1 . 

[0156] If the response is negative, that means that all the locations between the 
positions Xj and x i+1 have not yet been considered. Step S64 is then followed by the 
previously described step S61. 

[0157] If the response is positive at step S64, that means that all the locations 
between the positions x i and x i+1 have been tested in order to know whether they 
have already been encoded or not. The length of vector V i has been reduced by the 
number of locations met that were already encoded. 

[0158] Steps S60 and S64 are followed by the step S65 which is the entropy 
encoding of the current vector V { . 

[0159] The encoding parameter AC xi+1 of the location x i+1 is set to the value 1 . 
[0160] The following step S66 is a test to determine whether the parameter i has 
the value n-1, i.e. whether all the vectors to be encoded have been processed. 
[0161] If the response is negative, then step S66 is followed by the step S67 at 
which the parameter i is incremented by one unit in order to consider a following 
vector. Step S67 is followed by the previously described step S55. 
[0162] When the response is positive at step S66, the encoding of the path 
according to the invention is terminated. 

[0163] Figure 9 represents the application of the vector reduction according to the 
first embodiment of the invention. 

[0164] The vector Vj determined at step S55 has a length of five locations between 
the locations x i and x i+1 . Three of the locations between x i and x i+1 have already 
been encoded and are represented by crosses. 

[0165] The result of the processing of vector Vj is to reduce the length of that 
vector by three, so as to skip the locations already encoded. The resulting vector 
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V'j has a length of two locations corresponding to the intermediate location not yet 
encoded and to the location of x i+1 . 

[0166] Figure 10 represents the application of the vector reduction to the vectors 
of Figure 8. 

[0167] The total amplitude of the vectors is 33 before processing according to the 
invention. It is 28 by virtue of the invention. 

[0168] The vectors are shorter, which reduces the information to be transmitted 
and which corresponds to an increase in compression. 

[0169] The second variant of step S5 of encoding the locations of the coefficients 
of the current block is detailed in Figure 1 1 in the form of an algorithm comprising 
steps S500 to S515. 

[0170] Step S500 is the determination of a path in the current block. 
[0171] In a similar manner to the first variant, a path is determined by an initial 
coefficient and the list of vectors joining the other coefficients. Each coefficient of 
the path different from the initial coefficient is represented by a vector describing 
its location with respect to the preceding coefficient in the path. It should be noted 
that the path does not necessarily include all the coefficients of the current block. 
This is because it is possible to encode only part of the coefficients and to set the 
other coefficients to the value zero at the time of later encoding. 
[0172] For example, Figure 12 represents a path comprising six coefficients of 
which the locations are connected by five vectors V x to V 5 . The locations of the 
coefficients of the path are located by their coordinates (x i9 y^, (x 2 , y 2 ), . . ., (x 6 , y 6 ) 
in a coordinate system of two dimensional space. 

[0173] In what follows, the locations {(x i9 y^, (x 2 , y 2 ) . . ., (x n , y n )} of coefficients 
in the current block are considered, where n is an integer equal to the number of 
coefficients of the path. 

[0174] The following step S501 is an initialization at which a binary encoding 
parameter AC xy is associated with each of the locations (x, y) of the current block. 
All the encoding parameters AC are set to the value zero at step S501. 
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[0175] For the following processing, the encoding parameter AC xy has the value 
zero if the location (x, y) has not yet been encoded, and has the value one once the 
location (x, y) has been encoded. 

[0176] The following step S502 is an encoding of the first location (x i9 y 4 ) of the 
path. The parameter AC xlyl , which is associated with it is set to the value one. 
[0177] The following step S503 is an initialization at which a parameter i is 
initialized to the value one. 

[0178] The following step S504 is a test to determine whether the coordinate x { is 
greater than the coordinate y { . 

[0179] If the response is negative, step S504 is followed by the step S505 which is 
a test to determine whether the coordinate x i+1 is greater than the coordinate Xj. 
[0180] If the response is positive at step S504, then step S504 is followed by the 
step S506 which is a test to determine whether the coordinate y i+1 , is greater than 
the coordinate y^ 

[0181] When the response is positive at step S505, that step is followed by the step 
S507 at which the current vector Vj is decomposed into two vectors X i and Y v 
[0182] The vector X { is equal to: x i+1 - x i and the vector Yj is equal to: y i+1 - y d . 
[0183] The following step S508 is the reduction of the vector Xj. That reduction is 
similar to the reduction of a monodimensional vector described in the first variant. 
The monodimensional vector which is considered here is between the points (x i9 
y i+1 ) and (x i+1 , y i+1 ). 

[0184] When the response is negative at step S505, or when the response is 
positive at step S506, that step is followed by the step S509 at which the current 
vector Vj is decomposed into two vectors X A and Y { . 

[0185] The vector Xj is equal to: x i+1 - x i and the vector Y i is equal to: y i+1 - y^ 
[0186] The following step S510 is the reduction of the vector Y v That reduction is 
similar to the reduction of a monodimensional vector described in the first variant. 
The monodimensional vector which is considered here is between the points (x i+1 , 

Yd and ( x i + i> y i+ i)- 

[0187] When the response is positive at step S506, that step is followed by the step 
S51 1 at which the current vector Vj is decomposed into two vectors X i and Y it 
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[0188] The vector X { is equal to: x i+1 - x d and the vector Yj is equal to: y i+1 - 
[0189] The following step S512 is the reduction of the vector That reduction is 
similar to the reduction of a monodimensional vector described in the first variant. 
The monodimensional vector which is considered here is between the points (x i9 
y i+1 ) and (x i+1 , y i+1 ). 

[0190] Steps S508, S510, and S512 are followed by the step S513 which is the 
entropy encoding of the current vector Vj. 

[0191] The encoding parameter AC xi+lyi+1 of the location (x i+1 , y i+1 ) is set to the 
value one. 

[0192] The following step S514 is a test to determine whether the parameter i has 
the value n-1, i.e. whether all the vectors to be encoded have been processed. 
[0193] If the response is negative, then step S514 is followed by the step S515 at 
which the parameter i is incremented by one unit in order to consider a following 
vector. Step S515 is followed by the previously described step S504. 
[0194] When the response is positive at step S514, the encoding of the path 
according to the invention is terminated. 

[0195] Figure 13 represents the application of the vector reduction according to 
the second variant of the first embodiment of the invention. 
[0196] The vector Vj determined at step S500 extends between the locations (x i9 
Yi) and (x i+1 , y i+1 ). The locations already encoded are represented by crosses. 
[0197] The vector Vj is decomposed into two vectors X i and Yj. In the case 
represented, it is the vector X { which is reduced (step S512). 

[0198] The length of the resulting vector Vj is reduced with respect to the original 
vector Vj. 

[0199] Figure 14 represents the application of the vector reduction to the vectors 
of Figure 12. 

[0200] The total amplitude of the vectors is 15.65 before processing according to 
the invention and 14.13 after processing. The amplitudes have been measured in 
Euclidean distance. 

[0201] The vectors are shorter, which reduces the information to be transmitted 
and which corresponds to an increase in compression. 
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[0202] Figure 15 represents a first embodiment of a method of decoding data 
previously encoded according to the method of Figure 3. 

[0203] This method is implemented in the decoding device and comprises the 
steps S20 to S27. 

[0204] The method is carried out in the form of an algorithm which can be stored 
in whole or in part in any means of information storage capable of cooperating with 
the microprocessor. This storage means is readable by a computer or by a 
microprocessor. The storage means is integrated or not into the device, and may be 
removable. For example, it may comprise a magnetic tape, a diskette or a 
CD-ROM (fixed memory compact disk). 

[0205] Step S20 is an initialization at which the encoding data of the first block 
are considered. 

[0206] The following step S21 is the decoding of the amplitude model. The result 
is the amplitude model of the coefficients of the current block. 
[0207] Step S20 is also followed by step S22 which is the decoding of the 
coordinates of the coefficients. This decoding will be detailed below. The path is 
decoded to provide the locations of the coefficients in the decoded block. 
[0208] Steps S21 and S22 are followed by step S23 which is the association of an 
amplitude with each decoded location. For this, a value provided by the amplitude 
model is attributed to each coefficient location determined at step S22. The 
location of the coefficient decoded receives the component of the code 
vector constituting the amplitude model. 

[0209] The following step S24 is the setting to zero of the locations of the block 
which were not in the path. The result of steps S23 and S24 is a block B'. 
[0210] The following step S25 is an inverse DCT transformation of the block B' 
obtained at the preceding step. 

[0211] The following step S26 is a test to determine whether the current block is 
the last block to decode. 

[0212] If the response is negative, this step is followed by the step S27 at which a 
following block is considered. Step S27 is followed by the previously described 
steps S21 and S27. 
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[0213] If the response is positive at step S26, then the decoding of the image is 
terminated. 

[0214] A first variant of the step S22 of decoding the path is detailed with 
reference to Figure 16 in the form of an algorithm comprising steps S220 to S236. 
This is the case in which the path is formed in a monodimensional vector. 
[0215] Step S220 is an initialization at which a binary decoding parameter AD X is 
associated with each of the locations x of the current monodimensional vector. All 
the decoding parameters ADx are set to the value zero at step S220. 
[0216] For the following processing, the decoding parameter ADx has the value 
zero if the location x has not yet been decoded, and has the value one once the 
location x has been decoded. 

[0217] The following step S221 is a decoding of the first location x x of the path. 
The parameter AD xl which is associated with it is set to the value one. 
[0218] The following step S222 is an initialization at which a parameter i is 
initialized to the value one. 

[0219] At the following step S223, the vector Vj is decoded. The location x i+1 = 
V x + x i is considered. A location parameter x is set to the value x^ 
[0220] The following step S224 is a test to determine if the current vector V x is 
oriented in one direction or the other. For example an axis of coordinates of which 
the origin is at x x is considered. Thus, as for the encoding, it is tested here whether 
the vector V i is positive or negative. 

[0221] When the current vector V { is positive, step S224 is followed by step S225 
at which the location parameter x is incremented by one unit, in order to consider 
the location immediately to the right of the current location. 
[0222] The following step S226 is a test to determine whether the decoding 
parameter of location x has the value one, that is to say whether location x has 
already been decoded. 

[0223] If the response is positive, step S226 is followed by step S227 at which the 
value of the location x i+1 is incremented by one unit, so as to "skip" the location 
already decoded. 



-24- 



[0224] If the response is negative at step S226, this step is followed by the step 
S228. In the same way, step S227 is followed by the step S228. 
[0225] Step S228 is a test to determine whether the current value of the location 
parameter x is equal to x i+1 . 

[0226] If the response is negative, that means that all the locations between the 
positions x { and x i+1 have not yet been considered. Step S228 is then followed by 
the step S225 already described. 

[0227] If the response is positive at step S228, that means that all the locations 
between the positions x i and x i+1 have been tested in order to know whether they 
have already been decoded or not. The location of the coefficient x i+1 has been 
deduced from the location of x { and from the reduced vector V i9 taking into account 
the locations encountered that have already been decoded. 

[0228] When the current vector Vj is negative, step S224 is followed by step S229 
at which the location parameter x is decremented by one unit, in order to consider 
the location immediately to the left of the current location. 
[0229] The following step S230 is a test to determine whether the decoding 
parameter of location x has the value one, that is to say whether location x has 
already been decoded. 

[0230] If the response is positive, step S230 is followed by step S231 at which the 
value of the location x i+1 , is decremented by one unit, so as to "skip" the location 
already decoded. 

[0231] If the response is negative at step S230, this step is followed by the step 
S232. In the same way, step S231 is followed by the step S232. 
[0232] Step S232 is a test to determine whether the current value of the location 
parameter x is equal to x i+1 . 

[0233] If the response is negative, that means that all the locations between the 
positions x i and x i+1 have not yet been considered. Step S232 is then followed by 
the step S229 already described. 

[0234] If the response is positive at step S232, that means that all the locations 
between the positions x { and x i+1 have been tested in order to know whether they 
have already been decoded or not. The location of the coefficient x i+1 has been 
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deduced from the location of x { and from the reduced vector V i? taking into account 
the locations encountered that have already been decoded. 

[0235] Steps S228 and S232 are followed by the step S233 at which the decoding 
parameter AD xi+1 of the location x i+1 is set to the value one. 

[0236] The following step S234 is a test to determine whether the parameter i has 
the value n-1, i.e. whether all the vectors to be decoded have been processed. 
[0237] If the response is negative, then step S234 is followed by the step E235 at 
which the parameter i is incremented by one unit in order to consider a following 
vector. Step S235 is followed by the step S223 already described. 
[0238] When the response is positive at step S234, that step is followed by the step 
S236 which is the transformation of the monodimensional vector into a 
bidimensional block. 

[0239] The second variant of step S22 of decoding the locations of the coefficients 
of the current block is detailed in Figure 17 in the form of an algorithm comprising 
steps S2200 to S221 1. This embodiment corresponds to the case of the encoding 
of vectors formed directly in the bidimensional code-block. 

[0240] Step S2200 is an initialization at which a binary decoding parameter AD xy 
is associated with each of the locations (x, y) of the current block. All the 
decoding parameters AD xy are set to the value zero at step S2200. 
[0241] For the following processing, the decoding parameter AD xy has the value 
zero if the location (x, y) has not yet been decoded, and has the value one once the 
location (x, y) has been decoded. 

[0242] The following step S2201 is a decoding of the first location (x i9 y^ of the 

path. The parameter AD xlyl which is associated with it is set to the value one. 

[0243] Furthermore, a parameter i is initialized to the value one. 

[0244] At the following step S2202, the vector V i is decoded. The vector V i is 

composed of a vector X { and a vector Y { . The location (x i+1 , y i+1 ) is considered such 

that: 

x i+1 = Xj + Xj. 

y i+ i = Yi + 
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[0245] The following step S2203 is a test to determine whether the coordinate x { is 
greater than the coordinate 

[0246] If the response is negative, step S2203 is then followed by the step S2204 
which is a test to determine whether the coordinate x i+1 is greater than the 
coordinate 

[0247] If the response is positive at step S2203, then step S2203 is followed by the 
step S2205 which is a test to determine whether the coordinate y i+1 is greater than 
the coordinate 

[0248] When the response is positive at step S2204, that step is followed by the 
step S2206 at which the vector X i is processed as a monodimensional vector (first 
variant). The monodimensional vector which is considered here is formed between 
the points (x i9 y i+1 ) and (x i+1 , y i+1 ). After increasing the monodimensional vector, 
the coordinate x i+i is then equal to: x { + X v 

[0249] When the response is negative at step S2204, or when the response is 
positive at step S2205, that step is followed by the step S2207 at which the vector 
Yj is processed as a monodimensional vector (first variant) The monodimensional 
vector which is considered here is formed between the points (x i+1 , y^ and (x i+1 , 
y i+1 ). After increasing the monodimensional vector, the coordinate y i+1 is then equal 
to: y. + Y^ 

[0250] When the response is negative at step S2205, that step is followed by the 
step S2208 at which the vector X i is processed as a monodimensional vector (first 
variant). The monodimensional vector which is considered here is formed between 
the points (x i9 y i+1 ) and (x i+1 , y i+1 ). After increasing the monodimensional vector, the 
coordinate x i+1 is then equal to: x i + X it 

[0251] Steps S2206, S2207 and S2208 are followed by the step S2209 at which 
the decoding parameter AD xi+lyi+1 of the location (x i+1 , y i+1 ) is set to the value one. 
[0252] The following step S2210 is a test to determine whether the parameter i has 
the value n-1, i.e. whether all the vectors to be decoded have been processed. 
[0253] If the response is negative, then step S2210 is followed by the step S221 1 
at which the parameter i is incremented by one unit in order to consider a following 
vector. Step S221 1 is followed by the previously described step S2202. 
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[0254] When the response is positive at step S2210, the decoding of the path 
according to the invention is terminated. 

[0255] Figures 18 to 21 concern the second embodiment of the invention. 
[0256] With reference to Figure 18, an embodiment of an encoding device 3a 
according to the second embodiment of the invention is adapted to encode a digital 
signal with the object of compressing it. The encoding device is integrated into an 
apparatus, which is for example a digital camera, a digital camcorder, a scanner, a 
printer, a photocopier, a fax machine, a database management system, or a 
computer. 

[0257] A digital camera la acquires a digital image IM. The image is transmitted 
to an encoding device 2a of which the operation will be detailed in what follows 
with the help of algorithms. 

[0258] The encoding device comprises a transformation module 21a. According 

to the invention, it comprises: 

means 22a for determining a number of samples to encode, 

means 23a for constructing a list comprising the determined number of 

samples, classified by decreasing amplitude. 

[0259] Finally, it comprises a module 24a for actual encoding. 

[0260] The encoded image may be transmitted by a transmission module 3a 

operating conventionally, to a decoding device 4a. As a variant, the encoded image 

is simply stored in memory to be decoded later. 

[0261] The decoded image IM' may for example be transmitted to a display device 
5a. 

[0262] Figure 19 represents a first variant of a method of encoding an image, 
according to the second embodiment of the invention. This method is 
implemented in the encoding device and comprises steps S la to SI la. 
[0263] The method generally comprises a transformation of the signal to be 
encoded, then the determination of an amplitude model of the coefficients output 
from the transformation. The locations of these coefficients are next encoded 
according to a method using a path established between the coefficients. 
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[0264] Such an encoding method is described for example in French patent 
application No. 01 06933, as follows. 

[0265] A first step will consist in an encoding of the sign of the coefficients. The 
coefficients being of real value, the sign information will first of all be encoded 
independently, by transmitting for example a binary value at 0 for a plus sign and 
at 1 for a minus sign. 

[0266] The rate is simple to measure since we have an entropic code which is 
associated with each displacement vector. The number of bits necessary for the 
encoding of a vector can thus be obtained directly. 

[0267] This type of problem is well known in combinatorics. It is known as the 
traveling salesman problem, in which a certain number of towns have to be gone 
through in an order which minimizes the total distance. In this case, our problem is 
a particular instance called the time-dependent traveling salesman problem. 
[0268] Indeed, the cost associated with the passage from one sample to another 
depends on the path already followed, since the decoding distortion of a sample 
depends on the rank at which it will be decoded. 

[0269] The solution of this type of problem is found by use of an evolutionary 
algorithm. Once the optimum path in terms of rate-distortion has been obtained, 
the first coefficient is encoded. 

[0270] It should be noted that it is necessary to encode the first coefficient 
transmitted in a different manner to the others. In order to simplify the encoding of 
this coefficient, one could for example arbitrarily encode its location using a vector 
of displacement with respect to the point of coordinates (16, 16), that is to say the 
center of the block. 

[0271] Next, for all the other coefficients except the last, the central processing 
unit will encode the displacement vector will be encoded. At this step, there are 
transmitted the coordinates of the vector of displacement from the coordinates of 
the previously encoded coefficient. 

[0272] Displacement vector is understood to mean the vector of which the 
coordinates are the differences between the coordinates of the second point and of 
the first. It should be noted that in theory this solution may seem less beneficial 
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than a conventional solution in which the position of the coefficients would be 
directly encoded since a difference over 32 coefficients is encoded over 6 bits 
instead of 5 (from -31 to +32). Nevertheless, the nature of the signal is such that in 
general, the coefficients of high amplitude are found to be accumulated in the same 
place in a sub-band. It is therefore worthwhile to encode the points with reference 
to each other, such that the displacement vectors remain very small on average. 
[0273] Thus, by applying an entropic encoding to these vectors, it will be possible 
to obtain a code that is more compact. 

[0274] For as long as all the vectors have not been encoded, the step of decoding 
the displacement vector for all the other coefficients except the last will be 
reiterated. Once all the coefficients have been encoded, the it is verified if all the 
blocks have been processed. 

[0275] The method is carried out in the form of an algorithm which can be stored 
in whole or in part in any means of information storage capable of cooperating with 
the microprocessor. This storage means is readable by a computer or by a 
microprocessor. The storage means is integrated or not into the device, and may be 
removable. For example, it may comprise a magnetic tape, a diskette or a 
CD-ROM (fixed memory compact disk). 

[0276] Step S la is a linear or non-linear transformation of a digital image IM to be 
processed according to the invention. 

[0277] In the preferred embodiment of the invention, the transformation is a 
discrete cosine transformation (DCT) by blocks, such as that applied under the 
JPEG standard. 

[0278] In variant form, another transformation is used, for example a discrete 
wavelet transformation, as in the JPEG2000 standard. 

[0279] The following processing is performed block by block, since the DCT 

transformation generates blocks. If the transformation does not generate blocks, 

the processing is applied globally to the entire image. 

[0280] The following step S2a is an initialization at which a first block is 

considered. 
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[0281] The following step S3a consists of classifying the coefficients of the 
current block according to decreasing amplitude. This results in a list P of 
coefficients. 

[0282] The following step S4a is the determination of an amplitude model. For 
this, a function for approximating the series of the classified coefficients is 
determined. This function is for example a decreasing exponential defined by a set 
of parameters which are determined by regression. French patent application No. 
01 06933 describes this step in detail, as follows. 

[0283] The method will choose, from the series of P coefficients arranged from the 
greatest to the least and named by the series x l9 x 2 , x 3 , . . . , x p , a set A comprising a 
finite number of R parameters a l9 a 2 , . . . , B R . For example, a decreasing 
exponential family could be chosen, written f A (x)=a 1 .exp(a 3 .x-i-a 4 )a 2 . 
[0284] The method will then determine among a family of real functions that 
which reproduces as accurately as possible the evolution of the values of the x { as a 
function of their rank i. 

[0285] At these steps, the parameters a 1? a 2 , a 3 , and a 4 will be determined such that 
the function f A is the closest possible to the points (1, x x ), (2, x 2 ), . . . ,(P, x p ). This 
typical problem is called function regression. Numerous techniques exist for 
finding the best set of parameters. For example minimization of the root mean 
square error RMSE could be sought, equal to 



[0286] To achieve this, it could be chosen to quantize each parameter a t into a 
finite number of possible values. There are then a finite number of different values 
of parameters of A. The RMSE value is then calculated for each set of parameters 
and is then stored with the value of the parameters associated. When all the 
parameters have been incremented, the set which minimizes RMSE from among 
the stored parameters will be chosen. 
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[0287] Figure 4 shows an example of an amplitude model A. A value A(k) 
supplied by the amplitude model corresponds to each integer value k along the 
x-axis. The value A(k) is an approximation of the amplitude of the k* coefficient 
classified in decreasing order. 

[0288] The following step S5a is the determining of an encoding cost associated 
with the list P. The cost of encoding a block is the function C = R + A.D, in which 
R represents the rate of transmission of the encoded form of the block, D represents 
the distortion generated in the block reconstructed after encoding and decoding, 
with respect to the original image, and A is a parameter of adjustment between 
compression of the image and distortion generated by the encoding. 
[0289] It should be noted that the distortion is calculated by performing the sum of 
the squared errors between the amplitude of each coefficient of the block and the 
amplitude which it will have after decoding. 

[0290] The following step S6a is the storage in memory of the list P as the 
optimum list P opt , if the encoding cost C calculated at the previous step is less than 
the encoding cost of the optimum list P opt previously stored in memory. When this 
step is gone through the first time, the current list is stored as the optimum list. 
[0291] The following step S7a is a test to determine whether the current list P has 
zero length. 

[0292] If the response is negative, this step is followed by the step S8a at which 
the current list P is modified. The last sample is withdrawn from the list, which 
results in a new list. Step S8a is followed by the previously described step S5a. 
[0293] If the response is positive at step S7a, that step is followed by the step S9a 
which is the encoding of the current block with the list P opt i.e. with the list which 
provides the lowest encoding cost. That step is also detailed in French patent 
application No. 01 06933, as follows. 

[0294] A first step will consist of an encoding of the sign of the coefficients. The 
coefficients being of real value, the sign information will first of all be encoded 
independently, by transmitting for example a binary value at 0 for a plus sign and 
at 1 for a minus sign. 
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[0295] Once the optimum path in terms of rate-distortion has been obtained, the 
first coefficient is encoded. 

[0296] It should be noted that it is necessary to encode the first coefficient 
transmitted in a different manner to the others. In order to simplify the encoding of 
this coefficient, one could for example arbitrarily encode its location using a vector 
of displacement with respect to the point of coordinates (16, 16), that is to say the 
center of the block. 

[0297] Next, for all the other coefficients except the last, the central processing 
unit will encode the displacement vector will be encoded. At this step, there are 
transmitted the coordinates of the vector of displacement from the coordinates of 
the previously encoded coefficient. 

[0298] Displacement vector is understood to mean the vector of which the 
coordinates are the differences between the coordinates of the second point and of 
the first. It should be noted that in theory this solution may seem less beneficial 
than a conventional solution in which the position of the coefficients would be 
directly encoded since a difference over 32 coefficients is encoded over 6 bits 
instead of 5 (from -31 to +32). Nevertheless, the nature of the signal is such that in 
general, the coefficients of high amplitude are found to be accumulated in the same 
place in a sub-band. It is therefore worthwhile to encode the points with reference 
to each other, such that the displacement vectors remain very small on average. 
[0299] Thus, by applying an entropic encoding to these vectors, it will be possible 
to obtain a code that is more compact. 

[0300] For as long as all the vectors have not been encoded, the steps of decoding 
the displacement vector for all the other coefficients except the last will be 
reiterated. Once all the coefficients have been encoded, it is verified if all the 
blocks have been processed. 

[0301] Step S9a comprises the encoding of the locations of the coefficients based 
on the list P t . To do this, a path is determined by an initial coefficient and the list 
of vectors joining the other coefficients. Each coefficient of the path different from 
the initial coefficient is represented by a vector describing its location with respect 
to the preceding coefficient in the path. It should be noted that the path does not 
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necessarily include all the coefficients of the current block. This is because it is 
possible to encode only part of the coefficients and to set the other coefficients to 
the value zero at the time of later encoding. 

[0302] Once the path has been determined, the coordinates of the initial coefficient 
are encoded by binary encoding and the vectors are encoded by entropic encoding. 
[0303] The encoded form of a block of the image comprises an amplitude model 
which supplies an approximation of the amplitude of the coefficients and a path 
which supplies an ordered series of the locations of the coefficients. The location 
of the k** 1 coefficient of this series is determined by the path and its amplitude is 
determined by the y-coordinate corresponding to the x-coordinate k according to 
the amplitude model. 

[0304] The following step S 10a is a test to determine whether the current block is 
the last block of the image to encode. 

[0305] If the response is negative, this step is followed by the step SI la at which a 
following block is considered. Step SI la is followed by the previously described 
step S3a. 

[0306] If the response is positive at step SlOa, then the encoding of the image is 
terminated. 

[0307] Figure 20 shows a second variant of a method of encoding an image, 
according to the second embodiment of the invention. This method is 
implemented in the encoding device and comprises the steps S20a to S28a. 
[0308] Here too, the method generally comprises a transformation of the signal to 
be encoded, then the determination of an amplitude model of the coefficients 
output from the transformation. The locations of these coefficients are next 
encoded according to a method using a path established between the coefficients. 
[0309] Such an encoding method is described for example in French patent 
application No. 01 06933, as follows. 

[0310] A first step will consist in an encoding of the sign of the coefficients. The 
coefficients being of real value, the sign information will first of all be encoded 
independently, by transmitting for example a binary value at 0 for a plus sign and 
at 1 for a minus sign. 
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[0311] The rate is simple to measure since we have an entropic code which is 
associated with each displacement vector. The number of bits necessary for the 
encoding of a vector can thus be obtained directly. 

[0312] This type of problem is well known in combinatorics. It is known as the 
traveling salesman problem, in which a certain number of towns have to be gone 
through in an order which minimizes the total distance. In this case, our problem is 
a particular instance called the time-dependent traveling salesman problem. 
[0313] Indeed, the cost associated with the passage from one sample to another 
depends on the path already followed, since the decoding distortion of a sample 
depends on the rank at which it will be decoded. 

[0314] The solution of this type of problem is found by use of an evolutionary 
algorithm. Once the optimum path in terms of rate-distortion has been obtained, 
the first coefficient is encoded. 

[0315] It should be noted that it is necessary to encode the first coefficient 
transmitted in a different manner to the others. In order to simplify the encoding of 
this coefficient, one could for example arbitrarily encode its location using a vector 
of displacement with respect to the point of coordinates (16, 16), that is to say the 
center of the block. 

[0316] Next, for all the other coefficients except the last, the central processing 
unit will encode the displacement vector will be encoded. At this step, there are 
transmitted the coordinates of the vector of displacement from the coordinates of 
the previously encoded coefficient. 

[0317] Displacement vector is understood to mean the vector of which the 
coordinates are the differences between the coordinates of the second point and of 
the first. It should be noted that in theory this solution may seem less beneficial 
than a conventional solution in which the position of the coefficients would be 
directly encoded since a difference over 32 coefficients is encoded over 6 bits 
instead of 5 (from -31 to +32). Nevertheless, the nature of the signal is such that in 
general, the coefficients of high amplitude are found to be accumulated in the same 
place in a sub-band. It is therefore worthwhile to encode the points with reference 
to each other, such that the displacement vectors remain very small on average. 
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[0318] Thus, by applying an entropic encoding to these vectors, it will be possible 
to obtain a code that is more compact. 

[0319] For as long as all the vectors have not been encoded, the step of decoding 
the displacement vector for all the other coefficients except the last will be 
reiterated. Once all the coefficients have been encoded, the it is verified if all the 
blocks have been processed. 

[0320] Here the invention concerns the initialization of the evolutionary 
calculation, and more particularly the determination of an initial population of lists 
of coefficients. This population comprises a predetermined number of individuals. 
[0321] The method is carried out in the form of an algorithm which can be stored 
in whole or in part in any means of information storage capable of cooperating with 
the microprocessor. This storage means is readable by a computer or by a 
microprocessor. The storage means is integrated or not into the device, and may be 
removable. For example, it may comprise a magnetic tape, a diskette or a 
CD-ROM (fixed memory compact disk). 

[0322] The encoding of a block of the image is considered here. The block 
considered comprises L coefficients. 

[0323] Step S20 is an initialization at which the number N of coefficients of the 
list to construct is initialized to the value L. 

[0324] The following step S21a is the construction of a list of coefficients which 
comprises the N largest coefficients of the block, classified in order of decreasing 
amplitude. 

[0325] The following step S22a is the addition of the current list P to the initial 
population of the evolutionary algorithm. 

[0326] The following step S23a is a test to determine whether the initial 
population comprises the desired number of individuals, which here are lists. 
[0327] If the response is negative, that step is followed by the step S24a which is a 
test to determine if the current list P is of zero length. 

[0328] If the response is negative at step S24a, that step is followed by step S25a 
which is a modification of the number N of coefficients of the list. For example, 
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the number N is reduced by 10 % and rounded to the nearest whole integer. Step 
S25a is followed by the previously described step S21a. 

[0329] If the response is positive at step S24a, it is no longer possible to construct 
lists with coefficients of the block in the manner set out above and that step is 
followed by the step S26a at which the population is completed with lists picked 
randomly. 

[0330] The following step S27a is the calculation of a list of coefficients by 
evolutionary calculation, based on lists constructed by the preceding steps as initial 
lists population. This calculation is performed as in French patent application No. 
0106933, as follows. 

[0331] The method will calculate the path which minimizes the cost of encoding, 
that is to say find the path linking Q elements of the set {x x , . . . ,x p } such that the 
cost of encoding of this path is minimum. A measurement of encoding cost is thus 
made available. 

[0332] In typical manner, this measurement is expressed, for a given signal, in the 
form B+AD, where B is the rate associated with the signal, D is the distortion with 
respect to the original after decoding of the signal, and A is a parameter fixed by 
the encoder or the user of the encoder. 

[0333] The parameter X determines the balance between compression and quality. 
The distortion is a simple measurement to calculate. This is because, when the 
decoder receives the information about the localization of the coefficient, the 
decoding value will be f A (k). 

[0334] Thus, at the encoder, when the k* sample is in course of being encoded, the 
distortion is measured between the value of this sample and f A (k). 
[0335] The rate is also simple to measure since we have an entropic code which is 
associated with each displacement vector. The number of bits necessary for the 
encoding of a vector can thus be obtained directly. 

[0336] This type of problem is well known in combinatorics. It is known as the 
traveling salesman problem, in which a certain number of towns have to be gone 
through in an order which minimizes the total distance. In this case, our problem is 
a particular instance called the time-dependent traveling salesman problem. 
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[0337] Indeed, the cost associated with the passage from one sample to another 
depends on the path already followed, since the decoding distortion of a sample 
depends on the rank at which it will be decoded. The solution of this type of 
problem is found by use of an evolutionary algorithm. 

[0338] If the response is positive at step S23a, this means that the number of lists 
of the initial population has been reached. That step is then followed by step S27a 
already described. 

[0339] Step S27a is followed by step S28a which is the encoding of the current 
block on the basis of the list obtained at step S27a. Here too, the encoding is 
performed as in French patent application No. 0106933, as follows. 
[0340] A first step will consist of an encoding of the sign of the coefficients. The 
coefficients being of real value, the sign information will first of all be encoded 
independently, by transmitting for example a binary value at 0 for a plus sign and 
at 1 for a minus sign. 

[0341] Once the optimum path in terms of rate-distortion has been obtained, the 
first coefficient is encoded. 

[0342] It should be noted that it is necessary to encode the first coefficient 
transmitted in a different manner to the others. In order to simplify the encoding of 
this coefficient, one could for example arbitrarily encode its location using a vector 
of displacement with respect to the point of coordinates (16, 16), that is to say the 
center of the block. 

[0343] Next, for all the other coefficients except the last, the central processing 
unit will encode the displacement vector will be encoded. At this step, there are 
transmitted the coordinates of the vector of displacement from the coordinates of 
the previously encoded coefficient. 

[0344] Displacement vector is understood to mean the vector of which the 
coordinates are the differences between the coordinates of the second point and of 
the first. It should be noted that in theoiy this solution may seem less beneficial 
than a conventional solution in which the position of the coefficients would be 
directly encoded since a difference over 32 coefficients is encoded over 6 bits 
instead of 5 (from -31 to +32). Nevertheless, the nature of the signal is such that in 
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general, the coefficients of high amplitude are found to be accumulated in the same 
place in a sub-band. It is therefore worthwhile to encode the points with reference 
to each other, such that the displacement vectors remain very small on average. 
[0345] Thus, by applying an entropic encoding to these vectors, it will be possible 
to obtain a code that is more compact. 

[0346] For as long as all the vectors have not been encoded, the steps of decoding 
the displacement vector for all the other coefficients except the last will be 
reiterated. Once all the coefficients have been encoded, it is verified if all the 
blocks have been processed. 

[0347] Figures 21 to 27 concern the third embodiment of the invention. 
[0348] With reference to Figure 21, an embodiment of an encoding device 3b 
according to the third embodiment of the invention is adapted to encode a digital 
signal with the object of compressing it. The encoding device is integrated into an 
apparatus, which is for example a digital camera, a digital camcorder, a scanner, a 
printer, a photocopier, a fax machine, a database management system, or a 
computer. 

[0349] A digital camera lb acquires a digital image EVL The image is transmitted 
to an encoding device 2b of which the operation will be detailed in what follows 
with the help of algorithms. 

[0350] The encoding device comprises means 22b for determining an amplitude 
model of the coefficients and a path between the coefficients. 
[0351] According to the invention, it further comprises prior means 21b for 
putting the coefficients in order as a function of their respective locations in the set 
of data. 

[0352] The encoded image may be transmitted by a transmission module 3b 
operating conventionally, to a decoding device 4b. As a variant, the encoded 
image is simply stored in memory to be decoded later. 

[0353] The decoding device comprises means 41b for decoding an amplitude 
model of the coefficients and a path between the coefficients, so as to form a first 
set of decoded coefficients. 
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[0354] According to the invention, it comprises means 42b for putting the decoded 
coefficients in order as a function of their respective locations in the first set of 
decoded coefficients so as to form a second decoded set. 

[0355] The decoded image IM' may for example be transmitted to a display device 
5b. 

[0356] The operation of the decoding device will be detailed below using 
algorithms. 

[0357] Figure 22 shows an embodiment of a method of encoding an image, 
according to the third embodiment of the invention. This method is implemented 
in the encoding device and comprises the steps Sib to S9b. 
[0358] The method generally comprises a transformation of the signal to be 
encoded, then the determination of an amplitude model of the coefficients output 
from the transformation. The locations of these coefficients are next encoded 
according to a method using a path established between the coefficients. 
[0359] Such an encoding method is described for example in French patent 
application No. 01 06933. 

[0360] The method is carried out in the form of an algorithm which can be stored 
in whole or in part in any means of information storage capable of cooperating with 
the microprocessor. This storage means is readable by a computer or by a 
microprocessor. The storage means is integrated or not into the device, and may be 
removable. For example, it may comprise a magnetic tape, a diskette or a 
CD-ROM (fixed memory compact disk). 

[0361] Step Sib is the division of the image to be processed into blocks. The 
blocks are adjacent and are for example of square form comprising 8x8 samples. 
[0362] The following step S2b is an initialization at which a first block is 
considered. 

[0363] The following step S3b is a linear or non-linear transformation of the 
current block. 

[0364] In the preferred embodiment of the invention, the transformation is a 
discrete cosine transformation (DCT) as for the JPEG standard. 
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[0365] In variant form, another transformation is used, for example a discrete 
wavelet transformation, as in the JPEG2000 standard. 

[0366] In the following step S4b, the coefficients of the current block B { are put in 
order as a function of their respective locations in the block. 
[0367] According to a first variant represented in Figure 23, a zig-zag path is 
followed through block Bj. 

[0368] In block B i the diagonals of coefficients oriented from bottom left to top 
right are considered. The path begins at the top left. The path passes through the 
diagonals either from bottom left to top right, or in the opposite sense. For a given 
diagonal, the direction of the path is the direction opposite to that of the preceding 
diagonal. 

[0369] Following this path direction causes the coefficients to be arranged in a 
one-dimensional vector V containing 64 coefficients. 

[0370] According to a second variant represented in Figure 24, step S4b comprises 
two sub-steps S41b and S42b. During sub-step S41b, block B; is divided into four 
blocks, termed sub-blocks, each being of 4x4 coefficients. 
[0371] The following sub-step S42b is the aggregation of the sub-blocks 
previously formed in a three-dimensional table. For this, the sub-blocks are 
considered in lexicographical order (from top left to bottom right). In this order, 
they are stacked up beginning at the highest level. The result is a 
three-dimensional table T with 4x4x4 coefficients. 

[0372] With reference again to Figure 22, step S4b is followed by step S5b which 
is a classification of the coefficients by decreasing amplitude. 
[0373] It should be noted that due to the DCT transformation, the coefficients of 
highest amplitude have more chance of being at the top left of the 8x8 block 
resulting from the transformation. 

[0374] Thus, in the first variant, the coefficients of greatest amplitude have more 
chance of being at the start of the vector V. The work of classification is thus 
partially made during the previous step, or at least is facilitated by that step. 
[0375] In the second variant, the coefficients of greatest amplitude have more 
chance of being at the highest level of table T. 
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[0376] The classification is performed in conventional manner, for example by 
bubble sorting. This results in a list P of coefficients classified by decreasing 
amplitude. 

[0377] The following step S6b is the determination of an amplitude model. For 
this, a function for approximating the list P of the classified coefficients is 
determined. This function is for example a decreasing exponential defined by a set 
of parameters which are determined by regression. French patent application No. 
01 06933 describes this step in detail, as follows. 

[0378] The method will choose, from the series of P coefficients arranged from the 
greatest to the least and named by the series x l9 x 2 , x 3 , . . . , x p , a set A comprising a 
finite number of R parameters a l9 a 2 , . . . , B R . For example, a decreasing 
exponential family could be chosen, written f A (x)=a 1 .exp(a 3 .x-i-a 4 )a 2 . 
[0379] The method will then determine among a family of real functions that 
which reproduces as accurately as possible the evolution of the values of the x { as a 
function of their rank i. 

[0380] At these steps, the parameters a 1? a 2 , a 3 , and a 4 will be determined such that 
the function f A is the closest possible to the points (1, x x ), (2, x 2 ), . . . ,(P, x p ). This 
typical problem is called function regression. Numerous techniques exist for 
finding the best set of parameters. For example minimization of the root mean 
square error RMSE could be sought, equal to 



[0381] To achieve this, it could be chosen to quantize each parameter a t into a 
finite number of possible values. There are then a finite number of different values 
of parameters of A. The RMSE value is then calculated for each set of parameters 
and is then stored with the value of the parameters associated. When all the 
parameters have been incremented, the set which minimizes RMSE from among 
the stored parameters will be chosen. 
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[0382] Figure 4 shows an example of an amplitude model A. A value A(k) 
supplied by the amplitude model corresponds to each integer value k along the 
x-axis. The value A(k) is an approximation of the amplitude of the k* coefficient 
classified in decreasing order. 

[0383] The following step S7b is the encoding of the locations of the coefficients 
in the vector V or in the table T. To do this, a path is determined by an initial 
coefficient and the list of vectors joining the other coefficients. Each coefficient of 
the path different from the initial coefficient is represented by a vector describing 
its location with respect to the preceding coefficient in the path. It should be noted 
that the path does not necessarily include all the coefficients of the current block. 
This is because it is possible to encode only part of the coefficients and to set the 
other coefficients to the value zero at the time of later encoding. 
[0384] It should be noted that in the first variant, the coordinates of the initial 
coefficient and vectors are one-dimensional, whereas in the second variant the 
coordinates are three-dimensional. 

[0385] Once the path has been determined, the coordinates of the initial coefficient 
are encoded by binary encoding and the vectors are encoded by entropy encoding. 
[0386] The encoded form of a block of the image comprises an amplitude model 
which supplies an approximation of the amplitude of the coefficients and a path 
which supplies an ordered series of the locations of the coefficients. The location 
of the coefficient of this series is determined by the path and its amplitude is 
determined by the y-coordinate corresponding to the x-coordinate k according to 
the amplitude model. 

[0387] The following step S 8b is a test to determine whether the current block is 
the last block of the image to encode. 

[0388] If the response is negative, this step is followed by the step S9b at which a 
following block is considered. Step S9b is followed by the previously described 
step S3b. 

[0389] If the response is positive at step S8b, then the encoding of the image is 
terminated. 
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[0390] Figure 25 represents an embodiment of a method of decoding data 
previously encoded according to the method of Figure 22. 
[0391] This method is implemented in the decoding device and comprises the 
steps S20b to S27b. 

[0392] The method is carried out in the form of an algorithm which can be stored 
in whole or in part in any means of information storage capable of cooperating with 
the microprocessor. This storage means is readable by a computer or by a 
microprocessor. The storage means is integrated or not into the device, and may be 
removable. For example, it may comprise a magnetic tape, a diskette or a 
CD-ROM (fixed memory compact disk). 

[0393] Step S20b is an initialization at which the encoding data of the first block 
are considered. 

[0394] The following step S2 lb is the decoding of the amplitude model. The 
result is the amplitude model of the coefficients of the current block. This 
decoding is carried out as set out in the French patent application No. 01 06933, as 
follows. 

[0395] The parameters of the approximation function associated with the block to 
be processed are read, and they are decoded if necessary. 

[0396] Once this operation has been accomplished, the method will attribute, to 
the coordinate coefficient obtained by the displacement vector, the value of the 
amplitude given by the position of the coefficient by virtue of the approximation 
function reconstituted from its parameters. The decoding of the sign information 
will also be carried out so as to be able to attribute its sign to each coefficient 
decoded. 

[0397] The following step S22b is the decoding of the coordinates of the 
coefficients. This decoding is carried out as set out in the French patent application 
No. 01 06933, as follows. The path is described in order to provide the locations 
of the coefficients in the vector V of the table T. 

[0398] The method will attribute, to the coordinate coefficient obtained by the 
displacement vector, the value of the amplitude given by the position of the 
coefficient by virtue of the approximation function reconstituted from its 
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parameters. The decoding of the sign information will also be carried out so as to 
be able to attribute its sign to each coefficient decoded. 
[0399] Displacement vector is understood to mean the vector of which the 
coordinates are the differences between the coordinates of the second point and of 
the first. It should be noted that in theory this solution may seem less beneficial 
than a conventional solution in which the position of the coefficients would be 
directly encoded since a difference over 32 coefficients is encoded over 6 bits 
instead of 5 (from -31 to +32). Nevertheless, the nature of the signal is such that in 
general, the coefficients of high amplitude are found to be accumulated in the same 
place in a sub-band. It is therefore worthwhile to encode the points with reference 
to each other, such that the displacement vectors remain very small on average. 
[0400] The following step S23b is the construction of the vector V of the table T. 
For each coefficient of which the coordinates have been decoded at the preceding 
step, the order of the coefficient in the path determines its amplitude, since the 
coefficient of the path receives the amplitude A(k) corresponding to the 
x-coordinate k according to the amplitude model. 
[0401] The result of step S23b is a decoded vector V or table T. 
[0402] The following step S24b performs the reverse operations to those of step 
S4b at the time of encoding. The coefficients of the vector V or of the table T are 
ordered as a function of their respective locations in the vector V or table T to form 
a decoded block BY 

[0403] According to the first variant represented in Figure 26, the coefficients of 
the vector V are considered one after the other, from first to last, and are arranged 
in the block following the same zig-zag path as for encoding. 
[0404] According to the second variant represented in Figure 27, the levels of the 
table T are considered starting with the highest level, and are then arranged in 
lexicographical order in the block B\. 

[0405] The following step S25b is an inverse DCT transformation of the block B' 1 
obtained at the preceding step. 

[0406] The following step S26b is a test to determine whether the current block is 
the last block to decode. 
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[0407] If the response is negative, this step is followed by the step S27b at which a 
following block is considered. Step S27b is followed by the previously described 
step S21b. 

[0408] If the response is positive at step S26b, then the decoding of the image is 
terminated. 

[0409] Of course, the present invention is in no way limited to the embodiments 
described and represented, but encompasses, on the contrary, any variant form 
within the capability of the person skilled in the art. 



